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EXECUTIVE  SUMMARY 


We  measured  the  vapor-phase  absorptivity  coefficient  of  />/.v-(2-chloroethyl) 
sulfide  (HD)  in  the  mid-infrared.  We  used  agent  filled  saturator  cells  suspended  in  a  temperature 
controlled  liquid  bath  to  generate  continuous  streams  of  HD  diluted  in  nitrogen.  The  cells 
were  sent  to  a  variable  path  White  cell  and  measured  using  a  high  resolution  research  grade 
Fourier  transform  infrared  spectrometer.  The  purity  of  the  feedstock  was  verified  by  gas 
chromatography  (GC)  and  NMR.  The  concentration  of  HD  in  the  vapor  was  determined  with 
a  gravimetric  method,  and  the  vapor  generated  was  sampled  periodically  by  GC.  Ten  spectra 
at  different  concentration-pathlength  products  were  processed  line  by  line  through  least 
squares  analysis  using  MatLab®  to  produce  the  absorptivity  coefficient  of  the  compound  and 
the  statistical  uncertainty  in  the  data.  Uncertainties  in  the  data,  expanded  to  a  confidence  interval 
of  2o  (P=0.95),  are  Type-A  (4%)  and  Type-B  (4.2%)  of  the  absorptivity  coefficient. 

In  this  report,  we  report  a  comparison  of  our  data  to  that  obtained  by  another  laboratory  using 
a  different  vapor  generation  method. 


3 


Blank 


4 


PREFACE 


The  work  described  in  this  report  was  performed  under  the  direction  of  the 
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additional  copies  from  the  Defense  Technical  Information  Center;  unregistered  users  should 
direct  such  requests  to  the  National  Technical  Information  Service. 
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VAPOR-PHASE  INFRARED  ABSORPTIVITY  COEFFICIENT 
OF  fl/S-(2-CHLOROETHYL)  SULFIDE 


1.  INTRODUCTION 

We  have  obtained  the  high  resolution  vapor-phase  absorptivity  coefficients  of  the 
vesicant,  />/s-(2-chloroethyl)  sulfide  (HD),  in  the  spectral  range  of  4000-550  cm  1  in  units  of 
(pmol/moiy'm  ',  and  additionally  made  an  initial  computation  of  the  uncertainties  of  the  data. 
Previous  efforts  at  the  U.S.  Army  Edgewood  Chemical  Biological  Center  (ECBC)  to  obtain  the 
absorptivity  coefficients  of  HD  date  back  to  the  1960s  and  were  performed  at  lower  resolution  on 
a  grating  spectrometer.1  Given  the  extensive  use  of  vapor-phase  infrared  spectra  for  the  standoff 
detection  of  chemical  warfare  agents  (CWAs),  the  need  for  current  data  is  apparent. 

For  HD,  the  Chemical  Abstracts  Service  Registry  Number  is  505-20-2,  the 
molecular  formula  is  C^CLS,  and  the  molecular  weight  is  159.08.  Synonyms  for  the 
compound  include  2,2’-dichloroethyl  sulfide,  l,l’-thio-6A-(2-chloroethane),  sulfur  mustard, 
Levinstein  mustard,  and  S-lost.  The  structure  of  HD  is  shown  in  Figure  1 . 


Cl 


Cl 


Figure  1 .  Structure  of  6/'.v-(2-chloroethyl)  Sulfide  (HD) 


fi/5-(2-chloroethyl)  sulfide  is  a  delayed  action  vesicant  CWA  that  was  first 
employed  in  combat  during  WWI.  It  has  a  freezing  point  of  14.45  °C  and  a  vapor  pressure  of 
0.069  Torr  (9.2  Pa)  at  20  °C.2 

2.  EXPERIMENTAL  SECTION 

2.1  Instrumental  Details. 

The  system  used  to  generate  the  continuous  vapor  stream  was  an  adaptation  of  the 
saturator  cell  method  developed  at  ECBC  for  measuring  the  volatility  of  CWA  related 
compounds.1 


The  method,  modified  to  generate  continuous  streams  of  chemical  compounds  for 
obtaining  quantitative  vapor-phase  infrared  spectra,  has  been  used  to  measure  the  absorptivity 
coefficients  of  benzene4  and  a  variety  of  CWA  related  compounds.'  References  4  and  5  describe 
the  experimental  setup,  as  well  as  the  data  collection  and  post-processing  in  more  detail.  The 
saturator  passes  a  stream  of  nitrogen  carrier  gas,  obtained  from  the  boil-off  of  a  bulk  liquid 
nitrogen  tank,  across  a  conical  alumina  wicking  mechanism  in  a  glass  holder  filled  with  the 
analyte.  A  saturated  vapor-liquid  equilibrium  of  the  analyte  on  the  downstream  side  of  the 
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saturator  cell  results  with  the  concentration  of  the  analyte  determined  by  the  temperature  of  the 
liquid  phase.  By  suspending  the  saturator  cell  in  a  constant  temperature  bath,  the  concentration 
of  the  analyte  can  be  predicted  by  its  vapor  pressure  at  the  temperature  of  the  bath.  The 
apparatus  used  in  the  Quantitative  Fourier  transform  infrared  (FTIR)  Laboratory  uses  a  Brooks 
Model  5850S  mass  flow  controller  to  maintain  a  constant  flow  to  the  saturator  cell,  along  with  a 
second  mass  flow  controller  to  add  diluent  to  the  stream,  providing  an  additional  means  of 
adjusting  the  concentration  of  the  compound  delivered  to  the  White  cell  of  the  FTIR.  Linearity 
of  the  S  series  mass  flow  controllers  is  adjusted  using  a  second  order  polynomial,  resulting  in 
accuracies  of  approximately  1%  or  a  better  rate  at  flows  >25%  of  full  scale. 

Spectra  were  obtained  with  a  Bruker  Model  IFS/66V  FTIR  spectrometer.  The 
instrument  was  equipped  with  both  deuterated  triglycine  sulfide  and  mercury-cadmium-telluride 
(HgCdTe)  detectors  and  is  capable  of  obtaining  spectra  with  a  maximum  spectral  resolution  of 
0.1 125  cm'1  (unapodized).  The  interferograms  were  recorded  from  15798-0  cm'1  with  a 
resolution  of  0.125  cm'1.  Absorbance  (log  base- 10)  spectra  were  processed  with  boxcar 
apodization  and  2X  zero  filled  to  obtain  a  data  spacing  of  0.0625  cm'1.  The  instrument  was 
equipped  with  a  variable  path  White  cell.  The  experimental  data  used  path  lengths  of  4.057  and 
8.024  m.  The  temperature  of  the  White  cell  was  maintained  at  23  ±  0.1  °C  through  the  use  of  a 
thermostatically  controlled  chamber  enclosing  the  spectrometer  and  cell.  Data  were  acquired  at  a 
speed  of  60  KHz  (HeNe  laser  zero  crossing  frequency)  using  the  HgCdTe  detector.  Single  beam 
spectra  of  the  CWA  were  ratioed  against  spectra  of  clean,  dry  nitrogen.  To  minimize  the  effects 
of  nonlinearity  in  the  detector,  the  interferograms  were  processed  using  the  proprietary  Opus® 
nonlinearity  correction  function.  All  interferograms  were  archived  enabling  further  post¬ 
processing  of  data. 

Temperature  and  pressure  data  were  recorded  using  the  National  Institute  of 
Standards  and  Technology  (NIST)  traceable  digital  manometers  and  thermometers,  and  all  data 
were  archived.  Concentration-pathlength  (CL)  products  were  computed  in  units  of  pmol/mol(m) 
(ppm-m).  A  differential  pressure  manometer  had  previously  been  used  to  measure  the  dynamic 
pressure  in  the  White  cell  with  gas  flowing  into  the  cell,  and  the  ambient  pressure  was  plotted 
versus  the  differential  pressure.  The  resulting  equation  was  used  to  correct  the  readings  from  the 
ambient  pressure  manometer  to  the  pressure  in  the  White  cell.  The  CL  data  were  corrected  to 
296  K  and  1.0132  X  10?  Pa  (760  Torr)  using  the  ideal  gas  law. 

2.2  Feedstock. 

The  material  used  to  generate  the  vapor  streams  for  the  experiments  was  obtained 
from  the  Chemical  Sciences  Team  at  ECBC,  with  a  reported  purity  of  99%  (NMR  and  GC-Mass 
Spectrometry). 


We  verified  the  purity  of  the  feedstock  using  GC  with  a  flame  ionization  detector 
(FID)  at  98.9%.  The  primary  impurity  was  1,4-dithiane  at  1.1%  by  weight.  We  additionally 
checked  the  purity  of  the  effluent  vapor  from  the  White  cell  by  thermal  desorption  GC-F1D.  In 
the  effluent,  1,4-Dithiane  was  observed. 
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3. 


RESULTS  AND  DISCUSSION 


Four  trials  were  run  to  obtain  spectra  at  10-CL  products.  A  trial  is  defined  as 
filling  and  weighing  the  saturator  cell,  suspending  it  in  the  bath,  applying  a  stream  of  nitrogen  for 
a  measured  time,  acquiring  several  spectra,  stopping  the  nitrogen  and  removing  it  from  the  bath, 
and  reweighing  the  saturator  cell  after  drying  the  exterior  surfaces  and  re-equilibrating  to  room 
temperature.  The  first  trial  was  conducted  with  the  saturator  cell  suspended  in  a  water  bath 
maintained  at  23.3  °C,  with  the  remaining  trials  at  24.0  °C. 

Trace  water  vapor  was  removed  from  the  spectra  by  spectral  subtraction. 

In  the  spectra,  1,4-dithiane  was  also  observed.  The  concentration  of  the  compound  in  the  vapor 
was  computed  using  the  absorptivity  coefficient  of  1,4-dithiane  obtained  in  our  laboratory,  and 
the  features  of  the  impurity  were  removed  by  spectral  subtraction.  Baseline  corrections  were 
performed  with  several  spectra  using  a  two  point  linear  subtraction.  The  correction  in  no  case 
exceeded  0.001  A  [/l=-logio(7)]. 

The  composite  spectrum  (absorptivity  coefficient)  was  computed  using  spectra 
with  CL  products  ranging  between  207-849  pmol/mol(m)  (corrected  to  296  K  and  101325  Pa). 

As  an  initial  check  of  the  quality  of  the  data,  Beer’s  Law  plots  of  two  spectral  lines,  1212.5  and 
720.6  cm'1,  were  calculated  using  MatLab®,  showing  that,  at  least  for  these  two  spectral  lines, 
the  data  appeared  to  be  well  fitted,  with  no  points  lying  outside  the  95%  confidence  limits  for 
either  a  repeated  set  or  a  repeated  single  X,  or  the  95%  confidence  limits  for  a  Grubbs  Test  for 
Outliers6  (Figures  2  and  3). 


1212.5  cm'1 


Figure  2.  Beer's  Law  Plot  of  1212.5  cm  1  Line  in  the  Vapor-Phase  Spectrum 
of  bis-(2-chloroethyl)  Sulfide 
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720.62  cm'1 


Figure  3.  Beer's  Law  Plot  of  720.6  cm'1  Line  in  the  Vapor-Phase  Spectrum 
of  6A-(2-chloroethyl)  Sulfide 


The  absorptivity  coefficient  (a)  and  uncertainty  (Type-A,  2a)  were  computed  line 
by  line  within  the  spectral  range  of  4000-550  cm'1  using  a  MatLab  program  written  in-house. 
Values  of  (A=-logT)  >  1 .5  are  normally  assigned  a  weight  of  zero.  Because  A  did  not  exceed 
0.615,  all  values  of  A  were  weighted  at  one.  Figure  4  is  the  plotted  absorptivity  coefficient  (a), 
and  Figure  5  is  the  Type-A  uncertainty  for  the  computed  spectral  range.  The  figures  are  plotted 
with  a  in  (pmol/molf'm1.  To  obtain  a  in  (mg/m2)'1,  multiply  the  values  in  the  mantissa  of 
Figure  4  by  0. 1 5 1 8.  This  factor  is  derived  from  eq  1  using  the  molecular  weight  of  HD  ( 1 59.08) 

m2  f  24.15 
mg  v  mw 


")  _  mol 
)  fMnol(m) 
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x  10'4 


Figure  4.  Absorptivity  Coefficient  of  6/s-(2-chloroethyl)  Sulfide 


x  10'5 


Figure  5.  Uncertainty  in  the  Absorptivity  Coefficient  of  />/s-(2-chloroethyl)  Sulfide  (Type-A,  2o) 
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Table  1  provides  the  absorptivity  coefficients  in  (pmol/molf'm  1  and  (mg/m2)'1 
for  selected  bands  in  units  of  wavenumber  and  micrometers  (pm). 


Table  1.  Absorptivity  Coefficient  of  6/s-(2-chloroethyl)  Sulfide  for  Selected  Bands 


Wavenumber,  cm'1 
(Wavelength,  pm) 

Absorptivity  coefficient,  (pmol/molf'm  1 
[(mg/m2)'1] 

1299.25 

1.141E-04 

(7.70) 

(2.19E-05) 

1212.41 

2.291E-04 

(8-25) 

(3.477E-05) 

1037.16 

3.460E-05 

(9.64) 

(5.25E-06) 

720.59 

7.1 17E-04 

(13.88) 

(1.08E-04) 

In  general,  expanded  type-A  uncertainties  were  2-3%  of  the  absorptivity 
coefficient,  as  seen  in  Figures  6  and  7.  Figure  6  is  a  plot  of  absorptivity  coefficients  (abscissa) 
and  fractional  uncertainty  (Type-A,  UA,  2a)  (mantissa).  Figure  7  is  a  plot  of  the  absorptivity 
coefficient  and  uncertainty  and  also  includes  a  best  fit  of  the  data  points  obtained  by  least 
squares,  which  is  an  approximation  of  UA.~  ax+b.  For  the  fitted  line  in  Figure  7,  the  coefficients 
are  a  =  1.72  X  10'2  and  b  =  2.03  X  10"6  (Table  2). 


Reprocessed 


Absorptivity  coefficient,  (p-mol/mol^m'1)  x  iq'4 


Figure  6.  Absorption  Coefficient  (Abscissa)  and  Type-A  Uncertainty,  Fractional,  2a,  for  HD 
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Uncertainty,  U  (Type  A,  2o)  (p.mol/mcil) 


Figure  7.  Absorption  Coefficient  (Abscissa)  and  Type- A  Uncertainty,  2a,  for  HD 
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Table  2.  Type-A  Statistical  Uncertainty  for  HD  Vapor-Phase  Absorptivity  Coefficient 


Type  A 

2a  ~  ma  +  b 

Slope 

Intercept 

m 

b 

1.72  X  10-2 

2.03  X  10-6 

Type-B  estimated  standard  errors,  their  sources,  and  the  combined  Type  A  and  B 
uncertainties  are  provided  in  Table  3.  The  expanded  combined  Type  B  uncertainty  was 
computed  using  eq  2 

Ab  =  (AL2+AT2+AP2+AFTIR2+ANL2+AMR2)i/2  (2) 

The  expanded  uncertainty,  Type-A  +  B,  was  computed  with  eq  3 

Ua+b=  (Ab2+Aa2)i/2*2  (3) 

The  sources  of  uncertainty  and  their  fractional  values  as  well  as  an  explanation  of 
the  symbols  in  eq  1  are  given  in  Table  4.  Among  the  Type  B  uncertainties,  the  purity  of  the 
vapor  dominates  at  0.0 1 5  ( 1  o). 


Table  3.  Uncertainties  in  Absorptivity  Coefficient  of  HD  from  ECBC  Data 
where  a  >  0.0000716  (pmol/moiy'm  1 


Symbol 

Fractional  deviation 

Source 

AL 

0.005 

Pathlength 

AT 

0.0006 

Temperature  of  White  cell 

AP 

0.0003 

Pressure 

AFTIR 

0.0005 

Drift  in  spectrometer 

ANL 

0.01 

Nonlinearity  in  detector 

AMR 

0.005 

Mass  rate 

AD 

0.0075 

Dilution  rate 

Apurity 

0.015 

Purity  of  feedstock 

Ab 

0.021 

Combined  type  B  (la) 

Aa 

0.020 

Type-A  deviation  (la) 

Ua+b 

0.057 

Combined  Type  A+B, 
expanded  to  2a 
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4. 


INTERLABORATORY  COMPARISON 


Comparison  of  data  between  laboratories,  especially  when  obtained  using 
different  methods,  can  be  useful  in  assessing  the  experimental  methods  and  for  validating  the 
accuracy  of  the  data.  One  other  source  for  the  high  resolution  absorptivity  coefficient  of 
/>A-(2-chloroethyl)  sulfide  for  comparison  to  the  data  acquired  at  ECBC  is  that  obtained  by 
Pacific  Northwest  National  Laboratory  (PNNL)  at  Dugway  Proving  Ground.  The  PNNL 
database  includes  a  detailed  discussion  of  experimental  methods  and  statistical  processing.  Each 
absorptivity  coefficient  spectrum  is  provided  with  a  metadata  file  listing  physical  data  for  the 
compound,  temperature  and  pressure  range,  number  of  individual  spectra  included  in  the  fit  of 
the  absorptivity  coefficient,  and  Type  A  and  Type  B  uncertainties. 

As  can  be  seen  in  Figure  8,  the  spectra  from  ECBC  and  PNNL  are  qualitatively 
similar.  Both  the  PNNL  and  our  laboratories  have  undertaken  an  extensive  effort  to  ensure  the 
wavelength  accuracy  of  the  Bruker  IFS-66V  spectrometers  used  to  acquire  spectra,  and  that 
deviations  average  <1  part-per-million  in  the  mid-infrared.  Indeed,  we  have  in  the  past  used 
spectra  of  compounds  exhibiting  significant  rotational  fine  structure  that  had  been  acquired  by 
PNNL  to  perform  spectral  subtraction  of  spectra  obtained  in  our  laboratory  with  few  residual 
features  remaining.  When  the  spectra  are  expanded  within  the  range  of  4000-550  cm'1,  it  is 
possible  to  see  differences  in  the  baselines  of  the  two  spectra  (Figure  9).  These  are  more 
prominent  in  the  vicinity  of  2400  cm  1  and  near  the  low  frequency  cutoff  at  600-550  cm  1  In  the 
fingerprint  region  (Figure  10),  there  are  peaks  in  the  PNNL  spectrum  from  a  possible  impurity  in 
the  vicinity  of  1000  cm1,  as  well  as  a  shift  in  the  frequency  of  the  CH2  scissoring  band  from 
1454  to  1445  cm'1. 
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Figure  8.  Absorptivity  Coefficient  of  Vapor-Phase  HD  from  ECBC  (red)  and  PNNL  (blue) 
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Figure  9.  Absorptivity  Coefficient  Spectra  of  HD  Vapor  from  ECBC  (red)  and 

PNNL  (blue).  Spectra  are  zoomed  to  show  the  differences  in  the  baselines, 
particularly  near  the  low  frequency  cutoff  and  -2400  cm'1. 
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Figure  10.  Absorptivity  Coefficient  Spectra  from  ECBC  (in  red)  and  PNNL  (in  blue). 

Spectra  are  zoomed  into  the  fingerprint  region.  Note  an  apparent  impurity 
at  -1020-920  cm  1  and  a  shift  in  the  frequency  of  the  CH2  scissoring  vibration 
from  1454  to  1445  cm1. 


One  prominent  impurity,  1,4-dithiane  (DTH),  present  in  both  the  PNNL 
composite  spectrum  and  the  raw  spectra  from  our  laboratory,  can  be  readily  identified  by  a 
strong,  narrow  band  at  2955.4  cm1.  We  recently  acquired  the  absorptivity  coefficient  spectrum 
of  1 ,4-dithiane  in  the  vapor-phase.  Before  computing  the  absorptivity  coefficient,  the  features  of 
the  DTH  were  removed  from  individual  spectra  of  HD  using  spectral  subtraction  and  adjusted 
CL  products  were  computed  for  the  individual  spectra  that  accounted  for  the  concentration  of  the 
impurity  in  the  spectra.  Correcting  the  PNNL  composite  spectrum  for  1,4-dithiane  reduced  the 
integrated  area  of  the  spectrum  by  1  %. 

Differences  in  the  baselines  of  the  spectra  from  ECBC  and  PNNL  may  arise  from 
the  7th  order  polynomial  baseline  correction  method  used  by  PNNL.  Because  we  have  rarely 
observed  baseline  changes  of  A  >  0.001-0.002  [A  =  -log(I/Io)],  we  have  taken  a  rather 
conservative  approach.  We  use  a  MatLab  baseline  correction  algorithm  employing  linear  least 
squares  to  fit  the  baseline,  generally  to  no  more  than  2  or  3  points.  Data  taken  in  our  laboratory 
show  this  to  be  satisfactory  even  in  cases  of  baseline  shifts  significantly  greater  than  we  typically 
observe.  After  performing  a  multipoint  linear  correction  on  the  PNNL  spectrum,  the  baselines  of 
the  two  spectra  were  sufficiently  similar  to  improve  the  qualitative  and  quantitative  comparison 
of  the  two  spectra,  although  differences  were  still  seen. 
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The  PNNL  infrared  database  includes  the  imaginary  refractive  index  [i(/0]  of 
liquid  HD  computed  from  the  quantitative  transmission  spectrum.  The  absorptivity  coefficient 
of  HD  can  be  obtained  using  eq  4 


a  - 


cm 1 

ln(10)(/«o/) 


(4  n{v)K{v) 


(4) 


The  absorptivity  coefficient  spectrum  thus  obtained  was  then  recomputed  to  units  of  (g/cm3)(pm) 
using  the  molecular  weight  of  HD,  1 59.2. 


When  the  PNNL  and  ECBC  vapor-phase  spectra  are  overlaid  (Figures  1 1  and  12), 
it  is  possible  to  observe  what  appears  to  be  a  correspondence  between  features  in  the  PNNL 
condensed  and  vapor-phase  spectra.  Differences  in  frequencies  of  absorption  bands  are  common 
between  vapor-phase  and  condensed-phase  spectra.  Typically,  more  electronegative  or  polar 
bands,  or  those  subject  to  hydrogen  bonding  are  shifted  to  the  blue  (i.e.,  higher  frequency,  lower 
wavelength)  in  the  vapor-phase,  with  smaller  shifts  in  less  polar  bands. 8,9  This  is  apparent  in  the 
spectra  of  HD,  in  which  the  C-Cl  stretch  is  observed  at  721  cm'1  in  the  vapor-phase  and  at 
702  cm'1  in  the  condensed-phase.  The  effect  on  the  less  polar  methylene  scissoring  vibration  is 
less  pronounced,  with  a  shift  of  only  12  wavenumbers.  The  PNNL  vapor-phase  spectra  were 
acquired  using  a  syringe  pump  to  introduce  the  liquid  compounds  into  a  heated  finger  swept  by 
carrier  gas  at  10-15  L/min,  potentially  partially  aerosolizing  the  materials. 

We  then  used  spectral  subtraction  to  remove  the  features  apparently  arising  from 
the  presence  of  condensed-phase  material  in  the  vapor-phase  spectrum  from  PNNL.  The  result 
of  the  subtraction  can  be  seen  in  Figures  13  and  14.  Although  the  result  is  not  ideal  because  of 
the  lower  resolution  of  the  liquid-phase  spectrum  vis-a-vis  the  vapor-phase  spectrum,  the  liquid- 
subtracted  PNNL  spectrum  is  qualitatively  more  similar  to  the  ECBC  spectrum.  Unassigned 
peaks,  however,  are  still  present  in  the  PNNL  spectrum  around  1000  cm'1  and  1390-1350  cm1. 
The  latter  vibrations  are  generally  assignable  to  methyl  (CH3)  deformations,  not  present  in  HD, 
which  has  only  CH2.  A  common  impurity  in  HD  is  2-chloroethyl  ethyl  sulfide  (CEES),  which 
has  a  terminal  methyl.  In  addition  to  a  band  at  -1390  cm1,  not  seen  in  the  vapor-phase  spectrum 
of  HD,  CEES  has  a  band  at  -975  cm  1  that  is  nearly  40  times  more  intense  than  seen  in  HD, 
which  may  account  for  some,  but  not  all,  of  the  residual  differences  in  the  PNNL  and  ECBC 
spectra. 
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Figure  1 1.  Vapor-Phase  Absorptivity  Coefficients  of  HD  from  PNNL  (green)  and 
ECBC  (red),  and  the  Liquid  Absorptivity  Coeffients  (blue).  The  C-CI 
stretch  is  seen  at  a  A  of  19  wavenumbers  in  the  vapor-phase  vis-a-vis  the 
condensed-phase. 


irbance/  Wavenumber  (cm-1)  Overlay  Z-Zoom  CURSOR 

Figure  12.  Vapor-Phase  Absorptivity  Coefficients  of  HD  from  PNNL  (red)  and 

ECBC  (light  blue),  and  the  Liquid  Absorptivity  Coefficients  (dark  blue). 
The  methylene  scissoring  vibration  is  seen  at  a  A  of  12  wavenumbers 
in  the  vapor-phase  vis-a-vis  the  condensed-phase. 
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Figure  13.  Vapor-Phase  Absorptivity  Coefficients  of  HD  from  PNNL  after 
Spectral  Subtraction  (green)  and  ECBC  (red),  and  the  Liquid 
Absorptivity  Coefficient  (blue). 
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Figure  14.  Vapor-Phase  Absorptivity  Coefficients  of  HD  from  PNNL  after 
Spectral  Subtraction  (green)  and  ECBC  (red),  and  the  Liquid 
Absorptivity  Coefficient  (blue).  The  peak  corresponding  to  CH2 
scissoring  is  now  seen  at  the  same  frequency  in  both  vapor-phase  spectra. 
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The  ratio  of  the  CH  stretch  (3000-2800  cm1)  is  proportionately  more  intense  in 
the  PNNL  spectrum  than  in  the  ECBC  spectrum.  Table  4  shows  a  comparison  between  the  ratios 
of  CH  and  C-Cl  bands  for  the  two  spectra.  The  instrument  used  in  the  PNNL  work  was  modified 
by  the  addition  of  a  second  aperture  in  order  to  correct  for  artifacts  arising  from  double 
modulation  and  heating  of  the  aperture  by  the  source,  which  can  lead  to  a  distorted  absorption 
line  shape  that  shows  a  “tailing”  to  lower  frequencies.10  The  authors  of  the  referenced  study 
noted  that  the  artifacts  generally  tended  to  increase  the  integrated  areas  of  peaks  in  the 
fingerprint  region  relative  to  those  at  higher  wavenumbers.  The  National  Institute  of  Standards 
and  Technology  has  used  an  unmodified  Bruker  IFS-66V  infrared  spectrometer  to  compile  a 
quantitative  database.1 1  An  intercomparison  of  data  for  12  compounds  between  NIST  and  PNNL 
showed  that,  with  one  exception,  deviation  in  the  integrated  band  strengths  was  <2.5%.  The 
authors  have  attributed  the  small  apparent  magnitude  of  the  artifacts  to  the  use  of  a  slightly  large 
aperture  used  for  acquiring  the  NIST  data  (2.0  vs.  1 .5  mm)  and  the  lower  throughput  induced  by 
the  use  of  a  multipass  gas  (White)  cell  in  the  NIST  study.  The  deviation  seen  in  Table  4  is,  in 
any  case,  much  larger  than  that  noted  for  the  unmodified  PNNL  instrument.  Furthermore,  the 
data  in  Reference  10  were  obtained  using  a  20-cm  cell,  and  those  in  Table  4  were  acquired  using 
a  multipass  cell.  Artifacts  arising  from  double  modulation  and  aperture  heating  thus  apparently 
cannot  account  for  the  deviations  noted. 


Table  4.  Ratios  of  CH  to  C-Cl  in  Vapor-Phase  Spectrum  from  ECBC  and 
Spectrally  Subtracted  Spectrum  from  PNNL 


PNNL 

ECBC 

CH 

0.01397 

0.0091  1 

C-Cl 

0.01924 

0.01728 

CH/C-C1 

0.726 

0.527 

The  band  at  1037  cm'1,  although  weak,  is  of  interest  because  of  its  proposed  use 
as  an  identification  and  quantitation  band  for  LIDAR  measurements.  This  band  is  present  as 
well  in  the  condensed-phase  spectrum  of  HD  and  is  accompanied  by  a  second  band  at  1021  cm 
about  half  the  intensity  of  the  higher  frequency  vibration,  which  can  also  be  seen  in  the  PNNL 
vapor-phase  spectrum  (Figure  15).  Spectral  subtraction  of  liquid  HD  reduced  the  intensity  of 
this  band  in  the  PNNL  spectrum  by  (0. 000040-0. 00003)/0. 000033  X  100  =  21%  (vis-a-vis 
0.000034  in  the  ECBC  spectrum). 

By  computing  the  integrated  areas  of  the  PNNL  vapor-phase  spectrum  before  and 
after  subtracting  the  liquid-phase  spectrum,  it  is  possible  to  estimate  the  proportion  of  liquid  in 
the  vapor-phase  spectrum.  The  integrated  area  of  the  uncorrected  spectrum  is  0.05 1 20;  the  area 
of  the  corrected  spectrum  is  0.0460:  0.05120-0.0460  =  0.0591,  corresponding  to  an  average 
pathlength  for  undiluted  HD  of  0.045  pm. 
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Figure  15.  Band  at  1037  cm'1  in  Vapor-Phase  Spectra  of  HD  from  PNNL  (red)  and 

ECBC  (light  blue)  and  Liquid-Phase  Spectrum  from  PNNL.  After  spectral 
subtraction,  the  intensity  of  this  band  is  reduced  by  1 8%. 


A  comparison  of  the  integrated  areas  in  the  ECBC  and  PNNL  spectra  before  and 
after  spectral  subtraction  of  1,4-dithiane  and  liquid  HD  is  shown  in  Table  5.  The  data  are  for  the 
range  4000-550  cm'1.  The  metadata  file  associated  with  the  PNNL  vapor-phase  absorptivity 
coefficient  of  HD  indicates  uncertainties  of  3.5%  (Type-A)  and  <10%  (Type-B).  The  difference 
calculated  for  the  integrated  areas  of  the  PNNL  (subtracted)  and  ECBC  data  seen  in  the  last 
column  of  Table  5  is  well  within  the  limits  of  the  combined  uncertainties  of  the  experimental 
methods 


Table  5.  Comparison  of  Integrated  Vapor-Phase  Absorptivity  Coefficient  Spectrum 
of  HD  from  PNNL  and  ECBC 


PNNL 

ECBC 

(PNNL- 
ECBC)/ECBC 
X  100 

PNNL 

(after  spectral 
subtraction) 

ECBC 

(PNNL- 
ECBC)/ECBC 
X  100 

0.05746 

0.04390 

30.9% 

0.0460 

0.04390 

4.8% 
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5. 


CONCLUSIONS 


The  vapor-phase  absorptivity  coefficient  of  6/s-(2-chloroethyl)  sulfide  (HD)  was 
measured  within  the  range  4000-550  cm1.  Combined  uncertainties,  expanded  to  2 o,  are  5.7%  of 
the  absorptivity  coefficient  at  intensities  >5%  of  the  most  intense  absorption  feature. 


27 


Blank 


28 


LITERATURE  CITED 


1 .  Barrett,  W.J.;  Dismukes,  E.B.  Infrared  Spectra!  Studies  of  Agents  and  Field 
Contaminants',  9299-2018-1 1;  Southern  Research  Institute:  Birmingham,  AL,  1969; 
UNCLASSIFIED  Report  (AD-395  378). 

2.  Buchanan,  J.H.;  Buettner,  L.C.;  Tevault,  D.E.  Ambient  Volatility >  of  Bis-(2- 
chloroethyl)  Sulfide',  ECBC-TR-580;  U.S.  Army  Edgewood  Chemical  Biological  Center: 
Aberdeen  Proving  Ground,  MD,  2007;  UNCLASSIFIED  Report  (AD-B328  723). 

3.  Tevault,  D.E.;  Keller,  J.;  Parsons,  J.  Vapor  Pressure  of  Dimethyl  methylphosphonate 
ECBC-SP-004;  U.S.  Army  Edgewood  Chemical  Biological  Center:  Aberdeen  Proving  Ground, 
MD,  1999;  UNCLASSIFIED  Report  (AD-E491  779). 

4.  Williams,  B.R.;  Ben-David,  A.;  Green,  N.;  Hulet,  M.S.;  Miles,  R.N.;  Samuels,  A.C. 
Validation  and  Support  of  a  Quantitative  Vapor-Phase  Infrared  Instrument  Facility’  and 
Generation  of  a  Library  of  Chemical  Warfare  and  Related  Materials  by  Fourier  Transform 
Infrared  Spectroscopy,  ECBC-CR-076;  U.S.  Army  Edgewood  Chemical  Biological  Center: 
Aberdeen  Proving  Ground,  MD,  2006;  UNCLASSIFIED  Report  (AD-A471  712). 

5.  Williams,  B.R.;  Samuels,  A.C.;  Miles,  R.W.;  Hulet,  M.S.;  Ben-David,  A.  ECBC 
Quantitative  Vapor-Phase  Infrared  Spectral  Database,  In  Proceedings  of  the  2007  Scientific 
Conference  on  Chemical  &  Biological  Defense  Research,  13-15  November  2007',  SOAR-07-20; 
Chemical,  Biological,  Radiological  &  Nuclear  Defense  Information  Analysis  Center: 
Gunpowder,  MD,  2008. 

6.  Grubbs,  F.  Procedures  for  Detecting  Outlying  Observations  in  Samples. 
Technometrics,  1 969,  /  /( 1 ),  pp  1-21. 

7.  DOE/PNNL  Infrared  Spectral  Library,  Release  1 1.0;  Department  of  Energy/Pacific 
Northwest  National  Laboratory:  Richland,  WA,  2006. 

8.  Griffiths,  P.R.;  de  Haseth,  J.A.  Fourier  Transform  Infrared  Spectrometry,  Wiley- 
lnterscience:  New  York,  2007. 

9.  Nyquist,  R.  A.  The  Interpretation  of  Vapor-Phase  Infrared  Spectra',  Sadtler  Research 
Laboratories:  Philadelphia,  1984. 

10.  Johnson,  T.J.;  Sams,  R.L.;  Blake,  T.A.;  Sharpe,  S.W.;  Chu,  P.M.  Removing 
Aperture-Induced  Artifacts  from  Fourier  Transform  Infrared  Intensity  Values.  Appl.  Opt.,  2002, 
41,  p  283 1 . 

1 1 .  Chu,  P.M.;  Guenther,  F.R.;  Rhoderick,  G.C.;  Lafferty,  W.J.  The  NIST  Quantitative 
Infrared  Database.  J.  Res.  Natl.  Inst.  Stand.  Techno/.,  1999.  104,  p  59. 


29 


